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Single impurity problem in iron-pnictide superconductors is investigated by solving 
Bogoliubov-de Gennes (BdG) equation in the five-orbital model, which enables us to distin- 
guish s-i — and S++ superconducting states. We construct a five-orbital model suitable to BdG 
analysis. This model reproduces the results of random phase approximation in the uniform 

case. Using this model, wc study the local density of states around a non-rnagnctic impurity 
and discuss the bound-state peak structure, which can be used for distinguishing s-| — and s++ 
states. A bound state with nearly zero-energy is found for the impurity potential I ~ 1.0 eV, 
while the bound state peaks stick to the gap edge in the unitary limit. Novel multiple peak 
structure originated from the multi-orbital nature of the iron pnictides is also found. 
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Much effort has been devoted to elucidate the struc- 
ture of superconducting gap functions in the recently dis- 
covered iron pnictides. However, there are still many 
debates on this issue. Experimentally, the measurement 
of the quasi-particle interference pattern^^ and the exis- 
tence of the half-intcgcr flux-quantum transition''^ sug- 
gest that there is an internal 7r-phase shift of the gap 
function in iron-pnictide superconductors. Furthermore, 
nodal behaviors observed in some families of iron pnic- 
tides^' indicate that the repulsive interaction plays im- 
portant roles and the pairing with 7r-phase shift of the 
gap function sounds reasonable. On the other hand, it 
is claimed that the robustness against impurity dop- 
ing cannot be consistent with the above mechanism.^' 
Many theoretical studies have been performed in order to 
explain the high temperature superconductivity in iron 
pnictides. Most of these thcioricis have concluded that 
there exists inter-band sign reversal of the gap function 
(s-|-_ state). Recently, high temperature superconductiv- 
ity without inter-band sign reversal (s++ state) is also 
suggested. Considering these situations, it is impor- 
tant and urgent to develop some theories which enable 

us to distinguish the s^ and s++ states. There have 

been several theoretical proposals for detecting the s+_ 
state. However, it is a difficult and subtle problem to dis- 
tinguish s^ and s++ state than to detect other uncon- 
ventional superconducting states such as d-wave pairing, 
since the symmetries of the gap functions are the same 
for s+_ and s_|_+ state. 

In this paper, we investigate a single impurity prob- 
lem since it can serve as a possible method to detect the 
sign change of the gap function. Actually, this problem 
has been discussed mainly in some simplified two band 
models for iron-pnictide superconductors. ^^~^^' However, 
as we show here, detecting the S4._ state is subtle prob- 
lem and the entangled nature of the multiple bands and 
the Fermi surface structure, which characterize iron pnic- 
tides, should be taken into account. In the first part of 
this paper, we construct a model which can be used in the 



real-space single impurity problem. Although the five- 
orbital Hubbard model has been studied quite often, it 
is not suitable for our purpose since it does not give a 
superconducting ground state in the mean field theory. 
Here we construct a model whose interaction terms are 
chosen so as to reproduce the results in the five-orbital 
Hubbard model calculated within random phase approx- 
imation (RPA) . The hopping integrals of this model are 
taken from the downfolded result. Then, the local density 
of states (LDOS) around a single non-magnetic impurity 
is calculated by solving Bogoliubov-de Gennes (BdG) 
equation for the obtained model. The results show clear 
formation of in-gap bound state, which turns out to be a 
good quantity for distinguishing s+_ and s++ states. We 
find that the dependence of the spectrum on the impu- 
rity potential strength is unique to the present system, 
and there are novel multiple peak structures in LDOS 
for a certain parameter range. 

Now, we explain the five-orbital model used in this 
study. Hamiltonian for the used model can be written as 

ij 17 ab a a 

i o aba'b' 

+ X/ X/ X/ d'aa'b'b^iaiT^lba^ja'^'^ib'a 
(ij) a- aba'b' 

+ X/ X/ XI sia'b'b^iaa'^jba'^ja'^'^ib'cT, (1) 
{{ij)) f aba'b' 

where r* represents the position of the impurity site, and 
{ij) and {{ij)) represent the nearest and next nearest 
neighbor pairs respectively. Indices a or 6 run through 
to 4, where to 4 correspond to d^^^-r'^, dzx, dyz, 
dx2-y2, and d^y orbitals in this order. Hopping integrals 
ta,b-i,j are same as those in the Table I of Kuroki et al.^^ 
For details of these hopping parameters, see Ref 8) and 
a brief description is also available in our previours pa- 
per. The effect of impurity is simply treated as a local 
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Fig. 1. Real space distribution of the order parameter A"^" "' (r) 
around the impurity site with / = 1.0 eV. (a) for the orbital 1 
(d^x) and (b) for the orbital 3 (rf^2_ 2).^'^) 



potential shift, /, in eq. (1). The orbital dependence of 
the potential shift and the long-range impurity effects 
are neglected here. 

Coupling constants g(*^s in eq. (1) are chosen so as 
to reproduce the superconducting gap functions in the 
multi-orbital Hubbard model calculated within RPA.^^^ 
It was shown that the gap function obtained in RPA can 
be well reproduced with short-range pairings up to the 
next nearest neighbor sites when it is written in the or- 
bital representation, instead of the band representation. 
This point is reflected in eq. (1), where it is written in the 
orbital representation and the interaction terms are kept 
up to the next nearest neighbor sites. It was also shown 
that the pairings of d^z^-r^ and dxy orbitals are less im- 
portant since the density of states (DOS) just around the 
Fermi energy mainly comes from d^x/yz and d^^-y^ ot- 
bitals.^^-' Furthermore, the inter-orbital pairings are less 
important than the intra-orbital ones in the s-wave chan- 
nel. Based on these features, we have reduced the number 
of necessary parameters. Using the effective interaction 
calculated within RPA {V^^^{k), not shown), we obtain 
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in the unit of eV. Note that g^^j^,,^ = gj^aab ^^'^ parameters 
not shown here are set to be zero, g^aaaa represents the 
onsite repulsion and gioLa (g^aaa) represents the effective 
attractive interaction between the nearest (next nearest) 
neighbor sites, mainly induced by the broad (tt, tt) (sharp 
near (7r,0)) peak in V^l^{k). 

Hamiltonian, eq. (1), is solved in the mean field ap- 
proximation, where only the BCS type decomposition 
is used neglecting the Hartree-Fock type contribution. 
Order parameters are assumed to be site-dependent and 
then we obtain the multi-orbital version of the BdG equa- 
tion. In actual calculations, we introduce the energy cut- 
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Fig. 2. Partial LDOS in the case of / = 1.0 eV on the several sites 
near the impurity, (a) shows partial LDOS at just on the impurity 
site, (b) on the nearest, (c) the next-nearest, and (d) the third- 
nearest-neighbor site, respectively. Lines in each figure represent 
pi{r,ui) + p2{r,uj) (sum of the contributions from d^x and dyz 
orbitals) and p3{r,uj) (contribution from d^2_y2 orbital). LDOS 
at a site far from the impurity is also plotted. 



off, Wc, necessary in the standard BCS theory, i.e., we 
assume that the pairing occurs only between the elec- 
trons with energies e satisfying |e — /i| < lOc, where /i is 
the chemical potential. Throughout this paper, we use 
oJc = 0.1 eV. This choice of ujc does not affect the fol- 
lowing discussion. We solve the BdG equation in the lat- 
tices up to 28x28 sites with an impurity at the center, 
and determine the site-dependent order parameters by 
iteration. The electron number is fixed at n = 6.1 by 
adjusting the chemical potantial, fi. In order to calculate 
physical quantities, such as LDOS, we use the "super- 
cell" method^e) to improve the numerical accuracy. In 
this "super-cell" method, the 28x28 lattice with an im- 
purity is treated as a "unit cell" and the whole system 
is composed of 13x13 repetition of this unit cell. The 
connection beween the unit cells is represented by the 
wave vectors k = (27rn2./13, 27rnj,/13). The LDOS in this 
"super-cell" method is given by 



p(r,a;) 



Pa{r,uj) = -- VlmG^„ (fe,( 

TT ' * ' 



(3) 



(4) 



where Pair, w) is the partial LDOS for the orbital a, and 
^ra.r'a'i^T^) '^^ retarded Green's function in real 
space for the lattice sites r and r' in the unit cell. In fol- 
lowing figures, the Dirac delta functions in the LDOS are 
replaced by the lorentzian function with the half width 
7 = 0.001 eV. 
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Fig. 3. LDOS for various impurity potential /. Lines in each figure represents LDOS at the site far from the impurity site, on the 
impurity site, nearest neighbor site, next nearest neighbor site, and 3rd nearest neighbor site. 



Figure 1 shows the obtained order parameter and its 
modulation around the impurity site. Here, we consider 
the case of / = 1.0 eV. A representative order parameter, 
A'„>„" " (r), defined as 

^n.n.n./ X V-/^ , - V ' V 0^^' (c^ ) 

'-^aa y I — l-^r' a,ra ^ / , B aa' a' a\^r' a' a'~^ra' ah 

r' r' a' 

(5) 

is plotted in Fig. 1. Here, means that the summation 
over the next nearest neighbor sites of r. The amplitude 
of A"^'^ "' (r) is suppressed near the impurity and we can 
estimate the coherence length to be about four or five lat- 
tice spacing. At a site far away from the impurity, where 
the bulk behaviors are expected, we have A^i^ /Agg'^- ~ 

-0.41, A5'i'^" /A§f- 0.30, A§!J"'*74Ag3"- ~ 0.54, 

and Agg^'^'/Agg"- ~ —0.50. Roughly speaking, these 
ratios are consistent with those obtained in RPA for 
n = 6.1 and Jh/U = 0.2.1^) Although the value of 
A^i^ ^ /A'gg"- is slightly overestimated compared with 
RPA, the present result captures most of the features 
of RPA result. A'^-^-{r), which is defined for the nearest- 
neighbor sites, shows a similar behavior as A";j"" (r), 
while A™*''*''(r) = Ara,ra shows an overshooting behav- 
ior, i.e., it exceeds the bulk value at some points around 
the impurity. 

Next, we show LDOS around the impurity for the 
case of / = 1.0 eV in Fig. 2. In each figure, the partial 
LDOS for dzx/yz ^nd d^2_j,2 orbitals are plotted. From 
Figs. 2(a)-2(c), we can see the clear formation of the 
impurity induced in-gap bound state around impurity. 
When we look at the obtained spectrum more closely, 
we find that the peaks appear as a pair, i.e., at for 
each orbital, and that the width of each pair, E, depends 
on the orbital. As a result, total LDOS, which is the sum 
of the partial LDOS, shows novel multiple peak struc- 



tures (see also Fig. 3(b)). This feature is characteristic 
in iron pnictides and can be captured only in our realis- 
tic five-orbital model. At the third nearest neighbor site 
(Fig. 2(d)), on the other hand, the in-gap bound state 
peaks become small and the bulk-like coherence peaks 
recover. 

Impurity potential dependence of the LDOS is summa- 
rized in Fig. 3 where the results for / ranging from —8.0 
eV to -1-8.0 eV are shown. In each figure, we plot LDOS 
at several sites around the impurity up to third nearest 
neighbor site as well as LDOS far away from the impu- 
rity. As a function of / (/ > 0), (Fig. 3(a-d)), we find that 
the bound state appears at the edge of the gap (/ = 0.5 
eV), and moves toward nearly zero energy (/ = 1.0 eV), 
and then, goes back to the gap edge (/ = 2.0 eV and 
/ = 8.0 eV). As a result, the bound state formation is 
most prominent at / = 1.0 eV. On the other hand, when 
/ is negative, (Fig. 3(e-h)), the bound-state energies stick 
to the gap edge and do not approach the zero energy. In 
this case, the amplitude of the peak just grows with in- 
creasing |/|. Comparing these cases, we can see that the 
impurity with negative / has much less effects than the 
impurity with positive / when |/| is relatively small. In 
contrast, when |/| is large and in the unitary limit, LDOS 
structures are similar for both positive and negative /'s. 
This is reasonable since the impurity with large |/| works 
as a site onto which electrons cannot hop, whichever the 
sign of / is. 

Experimentally, LDOS can be measured by scanning 
tunneling microscopy/spectroscopy (STM/STS). Since 
the in-gap bound state around a non-magnetic impurity 
does not appear in s++ state, the bound state formation 

discussed above can be used to distinguish s^ and s_|_^ 

state. In particular, if the impurity potential is around 
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Fig. 4. Real space map of LDOS for (a) / = 1.0, ui = -0.002 eV 
and (b) / = 8.0, uj = -0.014 eV. 



J ~ 1.0 eV, the difference between s+_ and s++ states 
becomes most prominent since the in-gap bound state 
appears at the near zero-energy and the peak is very 
large. Even if the large in-gap state is not observed, the 
spectrmxi at the gap edge can be carefully investigated 

to distinguish the s.) and s++ state. Up to now, the 

superconducting gap has been successfully observed in 
Fe(TeSe) system,^' ^''^ and the detailed comparison of the 
theory and experiments is an interesting future work. 

We plot the real space map of LDOS in Fig. 4. In order 
to see the typical behavior of the impurity induced bound 
states, we choose the two cases with (a) / = 1.0 eV, ui = 
-0.002 eV and (b) / = 8.0 eV, lo = -0.014 eV, in which 
the bound state peaks in LDOS are clearly seen. For both 
cases, we find that the real space distribution of LDOS or 
the extension of the bound state is nearly isotropic. This 
will be due to the following two reasons: 1) Fermi surfaces 
of the used model are circular, and 2) the obtained gap 
function is nodeless in the bulk limit, which means that 
there is no special direction associated with the node. 
This feature will be changed if we consider a parameter 
region where the nodal gap function is obtained, or if the 
square Fermi surfaces appear. 

Here, we discuss the obtained results in connection 
to the T-matrix approximation, which has been used to 
study the formation of the bound state in the supercon- 
ducting state. ^^'^^^ This approximation applied on the 
similar five orbital model gives the results consistent with 
the present paper, although the site-resolved informa- 
tion is difficult to obtain in the T-matrix approxima- 
tion. We discuss the following two points. One is about 
the / dependence of the spectrum, which is understood 
from the particle-hole asymmetry of DOS in the present 
model. To be more specific, the DOS has larger weight 
below the Fermi energy than above the Fermi energy. We 
can show that this asymmetry results in the formation 
of the nearly zero-energy bound state when / is small 
and positive. (Details will be published elsewhere.) For 
other values of /, we can show similarly that the bound 
states are formed at the gap edge. The second point is 
about the multiple peak structure found in the case of 
/ = 1.0 eV. In this case, we can show that the T-matrix 
can be block diagonalized with each block having only 
intra-orbital contribution. This means that each orbital 
can sustain the bound state independently, and this is 
the origin of the multiple peak structure. 



It is tempting to speculate that the present results 
are connected to the experimentally observed robustness 
against the impurity doping. Namely, the impurity effect 
is weak when / is negative and small, which may explain 
the robustness against the impurity doping, although the 
magnitude of / and the relation to the residual resistivity 
should be taken account of carefully. 

In summary, we have calculated the LDOS around a 
non-magnetic impurity in the effective model for iron- 
pnictide superconductors. The model used here has re- 
alistic band structure and the effective interaction is de- 
termined so as to reproduce the RPA results of the gap 
function. It has been shown that the observation of the 
in-gap bound state enables us to distinguish the s++ 

and s_| states. An impurity with / ^ 1.0 eV gives 

large low energy peak, and we have also shown that the 
bound state peak appears at just below (or above) the 
gap edge in the unitary limit. An impurity with nega- 
tive and small / has relatively small effects on LDOS. 
The modeling of the case having the nodal gap function 
and of the recently proposed phonon-assisted orbital- 
fluctuation-medeated superconductivity are interesting 
extensions of the present study. Similar analysis on the 
quasi-particle interference patterns would also give im- 
portant information. 
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